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Abstract 

We study the equation of state of hyperon-rich matter for neutron stars using 
an extended relativistic mean-field model. We take special care of the recently 
proposed non-linear behaviour of the vector field which gives a much better 
description of Dirac-Briickner calculations. The hyperon-hyperon interaction 
is also implemented by introducing additional meson exchanges. These new 
terms avoid the instability found at high densitites in previous works while 
keeping the excellent description for finite nuclear systems. We also demon- 
strate within the mean-field approach that the presence of hyperons inside 
neutron stars on one hand and the hyperon-hyperon interactions on the other 
hand make the onset of kaon condensation less favourable. 

I. INTRODUCTION 

Nuclear matter at high densities exhibits a new degree of freedom: strangeness, hyperons 
and possibly kaons, occur at a moderate density of about 2-3 times normal nuclear density 
in neutron star matter []T[. These new particles influence the properties of the equation of 
state of matter and the global properties of neutron stars. Usually the relativistic mean 
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field (RMF) model is adopted for describing the equation of state with hyperons. But the 
standard approach suffers of several shortcomings: First, the parametrizations adjusted to 
reproduce the properties of nuclei fail at high densities due to some instabilities of the scalar 
selfinteraction (a possible solution to this problem has been suggested by Reinhard 0). 
Second the equation of state is much stiffer compared to nonrelativistic or Dirac-Briickner 
calculations and the effective mass is much smaller (m*/m « 0.5 — 0.6) in the ground state. 
A possible way out is to introduce a quartic selfinteraction term for the vector field |§. 
This model gives a good description of the properties of nuclei f|] and is in reasonable 
agreement with Dirac-Briickner calculations ||. In addition the instability found in the 
standard approach disappears. The equation of state for neutron stars (without hyperons) 
is considerably softened [|J. Third the hyperon-hyperon interaction becomes important for 
hyperon-rich matter present in the dense interior of neutron stars (note that nearly equal 
amounts of hyperons and nucleons are predicted jl]). The standard RMF model is not 
suited to reproduce the strongly attractive hyperon-hyperon interaction seen in double A- 
hypernuclei. An improved Lagrangian incorporating an additional pair of (hidden) strange 
meson fields remedies the situation |B|,[7J . These additional interaction terms have never been 
applied for the equation of state of neutron star matter before. 

In this paper we study the equation of state of neutron star matter for representative 
parameter sets which cover more or less all presently availabe fits of the relativistic mean 
field model to the properties of finite nuclei. The model is extended to include hyperons in 
a controllable way by fitting the parameters to hypernuclear data. In addition the hyperon- 
hyperon interaction is introduced through (hidden) strange meson exchanges. 

When describing the properties of dense matter it is also necessary to study the possibility 
of pseudoscalar meson condensation. In particular, much attention has been paid in recent 
years to the kaon condensation in neutron stars [Q. Most recent calculations based on 
chiral perturbation theory P,|l0"| show that kaon condensation may set in at densities only 
3-4 times larger than normal nuclear density. Nevertheless, these calculations do not take 
into account the presence of hyperons which will already occupy a large fraction of matter 
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when the kaons possibly start to condense. On the other side, the calculations including 
hyperons have not taken into account the possible kaon condensed phase. Only recently, 
some work has been started to incorporate the hyperon and kaon degrees of freedom at the 
same time pi|JI^| . In ref. [T1J the chiral Lagrangian by Kaplan and Nelson || is used for 



the baryons and kaons. The hyperons have not been included explicitly as constituents of 
the ground state, they are only considered through particle-hole excitations induced by the 
p-wave kaon-hyperon interaction. In this approach kaon condensation is predicted around 
3po- The other work [|l2j uses the standard relativistic mean field (RMF) approach for the 



baryon sector and the kaon-baryon interactions from the Kaplan-Nelson Lagrangian. It is 
shown that the critical density for kaon condensation is shifted to higher densities when 
hyperons are included (p c > 4p )- It has been observed that the nucleon effective mass is 
shifted towards negative values at high densities due to the kaon condensate. The authors 
conclude that the RMF model breaks down at these densities and that one has to go beyond 
the mean-field level. Thus, two approaches using the same kaon-baryon interactions come 
to different conclusions. In view of the present interest we also include the kaons in our 
RMF model by using kaon-baryon interactions motivated by meson exchange models. We 
demonstrate that the additional (hidden) strange meson fields mentioned above make kaon 
condensation less favourable even at very high densities. 

II. THE EXTENDED RMF MODEL 

We start from the standard Lagrangian 

B Z 

- gaB^B^BV ~ E 9uB^Bl^BV^ - £ Q p B^> ' Bl^B^ B R p (1) 
B B B 

where the sum runs over all baryons of the baryon octet (p,n,A,£ + ,£ ,E _ ,S°,2~). The term 
U (a) stands for the scalar selfinteraction 



U(a) = \mW + \f + Y ( 2 ) 



introduced by Boguta [O] to get a correct compressibility of normal nuclear matter. The 



parameters of this Lagrangian have been fitted to the properties of finite nuclei [fl^|I5]. It 



turned out that the best fits are obtained for the parameter sets with c < 0. In this case 
the functional form (Q) of the scalar field potential leads to an instability at high densities. 
Its traces are already seen in the nucleus 12 C 0. Hence, another stabilized functional form 
has been given by Reinhard || which eliminates the instability while keeping the good 
description of nuclei, especially the spin-orbit splitting. Alternatively, Bodmer proposed an 
additional selfinteraction term for the vector field of the form 

C v , = ±d(V,Vn 2 ■ (3) 

This modification leads to a nice agreement with Dirac-Briickner calculations at high den- 
sities ||. The reason is that the vector field is then proportional to p 1 / 3 in contrast to 
the linear dependence in the standard model. The fits to the properties of nuclei are quite 
succesful and the instability due to the scalar selfinteraction vanishes @]. 

The implementation of hyperons is straightforward. The corresponding new coupling 
constants have been fitted to hypernuclear properties |Uj. It turns out that the two coupling 
constants of the A (g aA and g u \) are strongly correlated because they are fixed by the depth 
of the A-potential 

U { A N) = g oA o*- + g^ (4) 



in saturated nuclear matter Hl7UT8|. Here denotes the potential depth of a baryon species 
i in matter of baryon species j. Hence one can use for example SU(6)-symmetry for the 
vector coupling constants 

g^JV = -3qjA = ~9ujS = 9ujE , 9 p n = -g P T, = 9pE , 9 P A = (5) 

and fix the scalar coupling constants to the potential depth of the corresponding hyperon. 
Following MM we choose 



U^ N) = U^ N) = -30 MeV , E4 JV) = -28MeV (6) 

in accordance with the available hypernuclear data. This constitutes our model 1. 

Nevertheless, these models are not able to reproduce the observed strongly attractive 
A A interaction irrespectively of the chosen vector coupling constant. The situation can be 
remedied by introducing two additional meson fields, the scalar meson / (975) (denoted as 
u* in the following) and the vector meson 0(1020) with the masses given in parenthesis ||[7J. 
The corresponding Lagrangian reads 

C YY = \ (d v a*d»a* - m>* 2 ) - + \m\^ 

- g^B^B^Bcr* - g^Bl^B^ ■ (7) 

B B 

The vector coupling constants to the 0-field are given by SU(6)-symmetry (see |7j for details) 

2^ 

2<70A — 2 fiV>£ = g<ps = — 9luN , g^N = . (8) 

The scalar coupling constants to the <r*-field are fixed by the condition 

UP w UP w 2Ui A) w 2U { A A) w -40 MeV (9) 

which is motivated by the one-boson exchange model D of the Nijmegen group and the 
measured strong AA interaction [0. Note that the nucleons are not coupled to these new 
fields. In the following we denote the extended model with hyperon-hyperon interactions as 
model 2. 

III. NEUTRON STAR MATTER WITH HYPERONS 

The equation of state for neutron star matter is derived by standard methods (see e.g. 
ref. [JTJ] for the standard RMF approach). The equations of motion for the meson fields in 
uniform matter are given by 

oa B 



2 * ( B ) 
m **° = 2^ 9a* BP S 
B 

mlV + dV 3 = J29^BP { v ] 



B 
B 
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B 



™|0O = E • (10) 



where ps and py denote the scalar and vector densities, respectively. The equations can be 
solved for a given total baryon density ps and charge density p c including the contributions 
from the free electrons and muons 



(B) 



Pv 

B 

pc = Y. ( ibpv ) + E w = o (ii) 

B l=e,fi 

where stands for the electric charge of a species i. In /^-equilibrium the chemical potentials 
of the particles are related to each other by 

Pi = Bi ■ p B + qi ■ p e (12) 

where Bi is the baryon number of a species i. This means that all reactions which conserve 
charge and baryon number are allowed, as e.g. 

n + n^A + n , A + A->H~+p , ... (13) 

Since we consider neutron stars in a long time scale, the strangeness quantum number is not 
constrained and the net strangeness is determined by the condition of /5-equilibrium. The 
above equations fixes already the fields and the composition of neutron star matter. The 
energy and pressure density can be derived from the grand canonical potential or the energy- 
momentum tensor in a standard way (see e.g. |J) and the generalization to the additional 
fields is straightforward: 

e = \mW + ^ + Y + \ m >* 2 + \< V * + l dV ° 

+ + \m\<& + £ « f F d^k^^f (14) 



-B 



6 



n 1 22 ^3 C 4 1 2 *2 , 1 2t/2,^t t/ 4 

r = — m^cr a a m„*a H — mr.Vn H — al/ n 

2 a 3 4 2 2 w 4 

+ + ±™w + e ^ r d3 *7ra (15) 



where 



m* B = m B + g aB cr + g a * B a* (16) 

is the effective mass of the baryon B (for leptons the on-shell mass is taken). 

In the following we study the composition of neutron star matter for the various param- 
eter sets obtained from fits to finite nuclei. First we take the sets NL-Z (which is set NL-1 



with a more consistent center-of-mass correction |I4]|), NL-SH |L5| for the Boguta- model and 
the sets PL-40 and PL-Z |2j with the stabilized scalar functional. The corresponding cou- 
pling constants are listed in table 1. The first parameter set shows a well known instability 
at high densities [0 which appears in neutron star matter around p ~ 0.5 fm -3 . At this 
point the effective mass of the nucleon gets negative due to the presence of the hyperons. 
For higher density no solution can be found. Negative effective nucleon masses appear also 
for set NL-SH at a density of p ~ 0.8 fm -3 for model 1 and at 1.0 fm~ 3 for model 2. The 
same holds for the case with the stabilized functional forms, sets PL-Z and PL-40 @, in 
model 1 and 2 around a similar density region. Note, that this happens even without the 



presence of a kaon condensate as found in [12]. One might wonder about the fact, that the 



early occurence of negative masses found here has not been seen in ||12|| . The explanation is 
simple: The standard parameter sets from fits to finite nuclei always favour a stiff equation 
of state to get a correct spin-orbit splitting which inhibits this instability as the fields are 
growing fast with the density. Fits to nuclear matter only can lead to a rather soft equation 
of state which avoids the instability due to a moderate rise of the fields. Nevertheless, the 
latter parametrizations are not able to describe the strong spin-orbit splitting seen in finite 



nuclei (for a detailed discussion see [Q) and therefore can not be considered as realistic. 
In the following we extended our calculation to densities beyond the instability, by taking 
always the absolut value of the effective mass (it is clear, that a much more refined 



procedure is necessary to treat the problem of negative effective masses, i.e. going beyond 
the mean- field appoximation). 

The behaviour for lower densities is quite similar in all these models. Fig. 1 shows as 
an example the composition of neutron star matter for the set PL-Z and model 2. The 
proton fraction raises rapidly and reaches maximum values over 20% around 2 — 3po- At 
this density, the hyperons, first A's and E~'s then S~'s, appear abundantly and the lepton 
fraction is considerably reduced. When the other hyperons are present at densities of 3 — 4po 
the number of A's even exceeds the number of neutrons, so that the dense interior resembles 
more a hyperon star than a neutron star [Q. Also the population of leptons gets negligibly 
low, as the electrochemical potential drops instead of growing. We have also checked for the 
occurence of Q~ which do not contribute in any of our calulations. 

The situation for higher densities differs for the various parameter sets mainly due to the 
instability mentioned above. Taking the absolute value of the effective nucleon mass still 
no solution can be found for the sets NL-Z and NL-SH above the critical density for model 
1. For the sets PL-Z and PL-40 we found a solution and the neutron star matter becomes 
pure hyperon matter dominated by A's and no nucleons appear above the critical density. 
Also the electrochemical potential changes sign and positrons and antimuons appear instead 
of electrons and muons. When turning on the hyperon-hyperon interactions (model 2) the 
hyperons feel at these high densities effectively an additional repulsion. Hence, nucleons 
are now present but still considerably less abundant than hyperons. Now a solution for set 
NL-SH can be also found. Nevertheless, these results should be taken with some care as we 
do not treat the problem of negative effective masses consistently. 

Now we turn our discussion to the case of models with vector selfinteractions. Two 
sets, named TM1 and TM2, exist so far in the literature 0] (see also table 1). It has 
already been found that the behaviour of the nuclear equation of state follows more nicely 
the trends of relativistic Briickner-Hartree-Fock calculations and that the maximum mass 
of a neutron star is reduced ||. Fig. 2 shows now the composition of neutron star matter 
for the set TM1 with hyperons and with hyperon-hyperon interactions (model 2). Up to 
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the maximum density considered, here 10po> ah effective masses remain positive and no 
instability occurs. The behaviour at moderate densities is quite close to that of the previous 
parametrizations: The proton fraction has a plateau at 2 — 4po and exceeds 20% which 
allows for the direct URCA process and a rapid cooling of the neutron star |I9| . Hyperons, 
first A's and £~'s, appear at 2p , then S~'s are populated already at 3p . The number of 
electrons and muons has a maximum here. The muons vanish when all the hyperons have 
been settled at p ~ 0.85 fm -3 while the electrons remain on the level of 2%. The fractions 
of all baryons show a tendency towards saturation, they reach asymptotically similar values 
of 8 — 15% which corresponds approximately to spin- and isospin-saturated matter. In the 
ideal case all baryons fractions would be the same due to the same spin-degeneracy factors. 

Switching off the hyperon-hyperon interactions, i.e. going from model 2 to model 1, we see 
again that negative effective masses occur at p = 1.3 fm~ 3 for set TM1 and at p = 0.95 fm -3 
for set TM2 due to the missing additional repulsive force for the hyperons. Nevertheless, 
the main population pattern is not changed considerably below p = 0.6 fm~ 3 when going 
from model 2 to model 1, but the number of hyperons, especially A's, are reduced at higher 
densities. Also the leptons vanishes for set TM1, but for the set TM2 positrons appear 
again, as the electrochemical potential changes sign above the occurrence of the instability. 
These results demonstrate the importance of the additional exchanges of a* and mesons. 
An additional repulsion is needed at higher densities to stabilize neutron star matter. 



The effective masses defined in (|16|) are plotted in fig. 3 for the various baryons for set 
TM1 and model 2. The effective mass of the nucleons reaches rapidly very small values 
and then saturates at higher densities never reaching negative values. For hyperons the 
behaviour is quite similar but less pronounced as their coupling constants to the cx-field are 
smaller than for nucleons. Nevertheless, the combined effect of the a and a* fields results 
in a rather constant gap between the effective masses over the whole density range shown. 
Note that the nucleons do not couple to the (hidden) strange scalar field a*. The overall 
small masses of the baryons at very high densities indicate that the neutron star matter 
approaches isospin-saturated matter with equal abundances, as can be seen in fig. 2. 



The pure scalar and vector field potentials are shown in fig. 4 (they have to be multiplied 
with the corresponding coupling constants given in table 1 to get the potential for a baryon 
species). In addition, we have also plotted the electrochemical potential. It reaches a 
maximum value of about 200 MeV around p « 2 — 3po and gets then consideraly lower at 
high densities instead of growing steadily. This behaviour is due to the onset of negatively 
charged hyperons (£~ and The vector potential exhibits initially a linear rise which 

then slows down due to the vector selfinteraction and the onset of hyperons at p ~ 0.3 fm~ 3 . 
At the same density the (hidden) strange fields are developed and reach quite high values for 
very dense matter. The 0-field shows a rather linear behaviour with density while the scalar 
fields seems to saturate at high densities. The isovector potential coming from p-exchange 
is quite small over the whole density range and never exceeds —60 MeV for the nucleons. 
The other potentials are one order of magnitude higher and have opposite signs. That also 
indicates a large cancelation between scalar and vector terms. The fields get so strong, that, 
for example, a 10% change of one coupling constant can alter the potential up to 100 MeV! 
Therefore, it is very important to fine tune the scalar coupling constants according to the 
potential depth (see eqs. (|B|) and (Q)). This has to be kept in mind for the forthcoming 
discussion of the properties of kaons in dense neutron star matter. 

The equation of state, the pressure density versus the energy density is given in fig. 5 for 
all parameter sets using model 2. For e < 800 MeV-fm -3 there exist mainly two bunches of 
curves. The parameter sets without vector selfinteraction are located in the upper branch. 
The sets TM1 and TM2 with vector selfinteraction are softer, i.e. they have a lower pressure 
for a given energy density and are therefore constituting the lower branch of curves. For 
higher energy densities, set TM2 shows to be stiffer than set TM1. For the upper branch one 
recognizes kinks in the curves which are due to the instability discussed above (i.e. negative 
effective masses). For model 1 (not shown), one sees a very pronounced jump in the curves 
due to the occurence of this instability. The equation of state gets then considerably softened 
reaching very low pressures. But for model 2 the jump is much less pronounced and even 
vanishes for the sets with vector selfinteraction. Another difference between model 1 and 
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model 2 is the slight softening of the equation of state for medium densities which is due to 
the attractive hyperon-hyperon interaction mediated by the <r*-meson exchange. For very 
high densities model 2 gets stiffer than model 1, because the contribution from the repulsive 
0-field overwhelms now the one coming from the attractive <r*-field which is saturating at 
very high densities. In addition the stiffest possible equation of state e = p is also drawn. 
The causality condition dp/de < 1 is fulfilled by all sets, so that the speed of sound remains 
lower than the velocity of light. The microscopic stability condition dp/de > is also 
satisfied except for the instability region. 

Finally, we want to study possible extensions of the RMF approach including other 
mesons. It is well known that pseudoscalar and pseudovector mesons do not contribute in 
the mean-field approximation. The same holds for all the kaon resonances as they couple off- 
diagonally. Tensor mesons also vanishes in infinite matter. Hence, the only remaining meson 
which can be added to the present Lagrangian is the scalar isovector meson ao(980) (the 
5-meson). To our knowledge this meson has been not considered so far for the equation of 
state of neutron stars. Fits to the properties of nuclei seems to indicate that its contribution 
are negligible for the discussion about binding energy, radius and surface thickness of stable 



nuclei [13]. Nevertheless, it might be important for very asymmetric systems. For the 
discussion in the next section it is quite useful to study the influence of this meson for the 
properties of neutron star matter. We implement the 5-meson in the standard way and 
choose the coupling constant from the Bonn-model (gsN — 5.95 p0[). The other coupling 
constants are scaled according to the isospin of the corresponding baryon (2g$N = g$z = 
2gs3 ? gsA = 0). As seen in fig. 4 the contribution from the vector isovector meson p is quite 
negligible compared to the other fields. Therefore, we expect only minor changes when 
introducing the scalar isovector meson 5. Indeed, the equation of state with the 5-meson 
gets only a little bit stiffer for lower density regions. For higher densities symmetric matter is 
reached also and the equation of state approaches the one without the 5-meson because the 
contribution of this meson simply vanishes in symmetric matter. The baryon composition 
shows some small changes for densities lower than p < 3po- As the 5-field is repulsive for 
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the protons, but attractive for the £~'s, the former ones are a little bit suppressed while 
the latter ones appear slightly earlier. These effects reduce the number of electrons and the 
maximum chemical potential is now about 150 MeV instead of 200 MeV. The meson fields 
do not change when the 5-field is introduced. The most pronounced effect, despite of the 
lower electrochemical potential, is the change in the effective masses of the baryons which 
are plotted in fig. 6. In contrast to fig. 3 every baryon has now a different effective mass. The 
biggest effect is seen for the S~ due to its large coupling constant to the 5-field. The effective 
mass of the S~ can even be lower than that of the A. The neutron has a significantly lower 
effective mass than that of the proton. For high densities it can be a relative difference of a 
factor of two! The absolute difference can be as high as 120 MeV. Nevertheless, the potential 
coming from the 5-field never exceeds 70 MeV for the nucleons and is therefore negligible 
compared to the rather high potential terms coming from the isoscalar fields (see fig. 4). 



IV. KAON INTERACTION IN THE RMF APPROACH 



The kaon-nucleon interactions has been studied recently within different approaches in- 



cluding the Bonn model pl |, the Nambu-Jona-Lasinio model [g2] and chiral perturbation 
theory [ffijfl . Antikaons in dense nuclear matter are of special interest because they can form 
a condensate f8[ leading to a considerable softening of the equation of state and a reduction 
of the maximum mass of neutron stars f9|. But this issue is still controversial. For in- 
stance rather strong nonlinear behaviour has been found by extrapolating KN-interactions 
derived from KN-scattering to dense nuclear matter in the NJL-model |22| where no hints 
for kaon condensation is seen. Applying the Bonn interaction scheme within the simple 
RMF approximation also leads a strong nonlinear density dependance disfavouring antikaon 
condensation pi| . Still, chiral perturbation theory gives a rather robust prediction for the 
onset of antikaon condensation at p c w 3 — 4po flOfl - 

We highly appreciate many efforts devoted to the consistent approach to kaon conden- 
sation based on the chiral perturbation theory pH|f2B] , |23| ] . But in our opinion this approach 
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has principal difficulties which may not be overcome in the near future. The main problem 
is that simple chiral models do not describe correctly the saturation properties of nuclear 
matter. By this reason the parameters used for the calculation of the kaon energy (coupling 
constants, effective masses, etc.) are not consistent with the properties of the background 
nuclear matter. Another and maybe related problem is that only linear terms in the baryon 
density are under control in the chiral approach. Therefore, all nonlinear terms are dropped 
out. In particular, the scalar density p s is identified with the baryon (vector) density p of 
the matter. This approximation certainly fails at high densities (for instance at p ~ 3po 
the difference can be about 80% [§]). Already using the scalar density, as it is dictated 
by Lorentz invariance, shifts the kaon condensation treshold to densities above 5po p4j , |26 



In ref. [p4|] it was also shown that many-body forces associated with the selfinteraction of 
the meson fields disfavour kaon condensation even more. By these reasons we use a more 
phenomenological approach based on the one-boson exchange model in the present paper. 

As it is well known, there exists an important difference between K + N and K~N- 
interaction. In contrast to K + N-scattering the situation with K~N- interactions is more 
complicated due to the existence of the A(1405)-resonance just below treshold. This makes 
the interpretation of the available K~N-scattering and K~-atomic data less transparent. 
Recently an improved fit of K~-atomic data was carried out assuming a nonlinear density 



dependence of the effective t-matrix \27\. It has been shown that the real part of the antikaon 



optical potential can be as attractive as 

C/| t « -200 ± 20 MeV (17) 

at normal nuclear matter density while being slightly repulsive at very low densities. This is 
in accordance with the low density theorem and K~p-scattering and effected by the A(1405)- 
resonance. Also another family of solutions have been found with a moderate potential depth 
around —50 MeV. Note that also the standard linear extrapolation gives only values of about 
—85 MeV [27]. These latter two solutions do not fulfill the low density theorem, i.e. getting 



repulsive at low densities. There exist some hints that the A(1405) is a quasi-bound state 
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in the t-channel pBLpTJ. It is not surprising then that the K N-scattering data can be 



explained by a simple vector meson exchange |ES] where the effects of this resonance come 



out automatically. This resonance seems to be less important in dense matter when the 



kaon mass is shifted down below mA(i405) — ~ 465 MeV. In ref. |29|] a separable potential 
was applied for the K~p-interaction for finite density. Indeed, it was found that the mass of 
the A(1405) is shifted upwards and exceeds the K~p threshold already at densities of about 
p ~ 0.4p - In this case the use of mean-field potentials may be justified. Most recently, the 
A (1405) has been also included in the chiral approach p5| . No effect has been found for 



the onset of kaon condensation |p5] , but note that the A(1405) has been only implemented 
as a new field contrary to the conclusion drawn in p8| , pl| . In a recent preprint the coupled 
channel analysis of Siegel and Weise p8|| has been applied for chiral perturbation theory 



]. The coupled channel formalism automatically generates the A(1405) and succesfully 
describes the low energy K~p-scattering data. 

In our previous work [^4j we have conjectured an important role which hyperons should 
play in determing the threshold of kaon condensation. Here we present the calculations 
which confirm these qualitative estimates. Recently, the effects of the presence of hyperons 



has been studied also in ref. [|T2J. The authors use the RMF model with scalar selfinteractions 
for the baryon sector and the chiral approach of Kaplan and Nelson |§ for the kaon-baryon 
interaction. It has been shown that the presence of hyperons, particular the S~, shifts the 
threshold for K~ condensation by several units of po- Note, however, that Ellis et al. do 
not take into account the terms which ensure that the S-term changes sign at low densities, 
as it should [[H]]. These additional terms will shift the onset of kaon condensation to even 
higher densities. 

In the following we adopt the meson-exchange picture for the KN-interaction simply 
because we use it also for parametrizing the baryon interactions. By analogy to the Bonn 



model |21| the simplest kaon-meson Lagrangian can be written as 

C K = d^Kd^K - m\KK - g aK m K KKa - g a * K m K KKa* 
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-guKKid^KVV-gpKKTKid^KW-g^Kid^K^ . (18) 

However, it is easy to see that this coupling scheme destroys the gauge invariance of the 
vector interactions represented by the equation 

3^ = . (19) 

From the equation of motion for the vector field one has 

dpF*" + mlV v = J2 9uB*b1 v ^b + 9uk{Ki d v K) (20) 

B 

and after applying a partial derivative on both sides 

m^V" = Y.^Bd^Bl^B) + g. K d v {Ki d v K) (21) 

B 

(the vector selfinteraction term is omitted here because it is purely phenomenological and 
the following arguments do not hold for it). The first term on the right hand side vanishes 
due to the baryon number conservation. The second term would also vanish for a free kaon, 
but it does not in the medium. Indeed, the conserved kaon current from the Lagrangian 
(H) reads 

■k ■ (f> dC dC \ 

= Kid^K - (id^K)K + 2KK{g u)K V lt + g pK rR p + g^) (22) 

so that this term does not vanish if the expectation value of the vector fields are nonzero. 
The situation can be corrected by introducing the additional terms 

C' K = C K + (g ujK V ll + g pK TR fl + g^ K( f) fl ) 2 KK . (23) 

The kaon-meson Lagrangian can now be written as 

C' K = D^KD^K - m 2 K KK - g aK m K KKa - g a * K m K KKa* (24) 

with the covariant derivative 
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D p = d p + ig uK Vp + ig pK rR p + ig^n ■ (25) 
The equation of motion for kaons in uniform matter then reads 

{<9 M <9 M + m 2 K + g aK am K + g a * K cr*m K + 2{g uK V + g pK r Ro,o + g^o)^ 

-{g^xVo + g pK T R 0t0 + g^x^o) 2 } K = . (26) 
The coupling constants to the vector mesons are chosen from the SU(3)-relations 

2g^ K = 2g pK = V2g<t>K = g^ P = 6.04 . (27) 

Decomposing the kaon field into plane waves one ontains the following dispersion relation 
for kaons (upper sign) and antikaons (lower sign) in uniform matter composed of nucleons 
and hyperons 

u K) r = \fm 2 K + m K (g aK a + g a * K o*) + k 2 ± (g wK Vo + g<p K (t>o + g P KRo,o) ■ (28) 

Note that due to the additional terms (^) the vector term appears linearly in the kaon 
energy. 

The coupling constant g a x could be taken from the Bonn model |2T| but here we adopt 
another prescription. We want to fix it to the value of the potential depth of kaons in the 
nuclear medium to study the strongly attractive potential seen in kaonic atoms [PTj] . The 
optical potential in symmetric nuclear matter (-Ro,o — 0, 0o = 0, no hyperons, k = 0) is in 
our approach given by comparison with eq. (p6j) 

2/i KN U* t = g aK om K - 2g u)K u R Vo - (g^xVo) 2 (29) 

where [Ikn is the reduced mass of the K~N-system. One notices that in general the optical 
potential and the energy shift (relative to the free mass) of the (anti)kaon do not coincide. 
In addition, the right hand side of eq. ( |29"D shows a nonlinear density dependance (see also 
|27||). At normal nuclear density the in-medium energy shift of the antikaon is about 20% 
less than the corresponding optical potential! In the following we study two possibilities: 
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we adjust the scalar coupling constant to U£L = — 120 MeV, as an upper bound for the 



standard fit and the lower family of solutions found in |p7jl , and to U^L = —200 MeV for 
the second family of solutions. 

For the KN case the low density theorem is not fulfilled due to the presence of the 
A(1405) resonance (in principle one has to use a coupled-channel formalism to get the correct 
scattering length which is beyond the scope of this paper). We argue here that the repulsive 
regime at very low density is of no consequence for the evaluation of K~ nuclear interaction 



at high densities as pointed out in p7|j29 



The onset of (s-wave) kaon condensation is determined by the condition 

- ii e = n K - = u K -(k = 0) . (30) 

When calculating fix- we have also taken into account the contribution from the (hidden) 
strange meson fields a* and 0o according to eq. (p8|). The coupling constant to the a* 
field has been taken from / -decay |31| (note that this value is obtaint despite of the fact 
that the mass of the fo is lower than the KK treshold which constitutes another source 
of uncertainty). While the <r*-field is attractive, the 0-field is repulsive for antikaons. The 
isovector field is also repulsive for the K~ as neutrons (and also E~) are more abundant than 
protons (S + ). As seen in fig. 5 the electrochemical potential decreases for higher densities. 
The relativistic effects, the reduced electrochemical potential, and the presence of the <fi- and 
the p-meson fields make the kaon condensation less favourable for strangeness-rich matter. 

In fig. 7 we plot ujk- — Hk- as a function of baryon density for model 2 and the various 
parameter sets considered. The upper curves correspond to U^£ t = —120 MeV whereas 
the lower curves to U^ t = —200 MeV. Below 2po the curves drop rather steep and linear. 
After the appearance of the hyperons, the slope of the curves changes dramatically and they 
even turn upward again. Antikaon condensation is only possible if the curves crosses zero. 
This does not happen at all for the upper curves and only occurs for set PL-40 (but in the 
unstable regime with negative effective masses!) and for set TM1 among the lower curves. 
Note that the two sets TM1 and TM2, which do not get negative effective masses and might 
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be the only two reliable calculations, give different predictions for the possible existence of 
a kaon condensed phase. While the former one reaches zero at p ~ 6po ; the latter one is at 
all densities at least 60 MeV above the critical value. 

For the case of U^ t = —200 MeV we have the problem that the optical potential of the 
kaon turns out to be attractive while experiment tells us that it is repulsive. Also the low 
density theorem is not fulfilled in our model. 

Hence, we discuss in the following the KN-scattering lengths, try to fix our parameters 
to the KN-scattering lengths and perform a G-parity transformation afterwards to get the 
KN-case. The isospin averaged scattering length in the tree approximation reads 

1 /=o . 3 I=1 m K I g a K9uN « 9u>k QuiN \ nn ,r ( / Q1 \ 

°o = 7 a o + 7 a o = ~a — n 7 V 9 2 9 — = -0-255 fm (31) 

4 u 4 u 47r (1 + m K m N ) \ w? a ml J 



where the experimental value has been taken from f52" |. It can be used to fix g a K for known 
9uK- Nevertheless, the presented Lagrangian is not able to get the correct scattering length 
for KN-scattering in the two isospin channels. A way out is, in addition to the p-meson, to 
introduce an additional coupling of the kaon to the scalar isovector meson which we will 
denote in the following as 5. The KN-scattering lengths for a given Isospin I on the tree 
level are then given by 

q /=1 _ m K / 9aK9aN g&KgSN _ ^ 9uK9wN _ ^9pK9pN_\ ^22") 

" An (1 + m^/mjv) \ m| ml m 2 p J 

1=0 m K ( 9aK9aN o gSKg5N n 9ujK9u)N . a 9pK9pN\ /qq n 

% u = - 7— / r 5 3 2 2 \-Q 9 • (33) 

47r (1 + m^jm^) \ m A a m & m^ m* J 

The importance of the 5-meson exchange contribution can be seen if one first sets it to zero 

and take SU(6)-symmetry for the vector coupling constants. Then the ag =0 scattering length 

depends only on the term coming from the a- meson exchange as g^Kg^N = ^g P Kg P N- This 

gives a large positive scattering length whereas experiment tells us that ag =0 = —0.09 fm is 



slightly negative |32|]. Taking g a ^ = 10, g^N = 13 as standard values for the RMF model (see 
table 1) one gets a 1 ^ w 0.4 fm without the 5-meson contribution. Including the 5- meson 
term and using g§N = 5.95 from the Bonn model [2D| one can fit both scatterings lengths 
nicely for 



g aK « 1.9 - 2.3 , g SK ~ 5.6 - 6.4 (34) 

for the various parameter sets used. A complete list of the coupling constants of the kaon 
can be found in table 2. Note that the values are surprisingly close to the one expected 
from the simple quark-model (g a K — 9un/3 ~ 3.3, gsx = gsN = 5.95). The effects of the 
5-meson on the properties of neutron star matter were discussed in the previous section and 
seemed to be quite insigificant. Hence, despite of the rather large coupling constant to the 
5-meson its influence on the kaon effective energy in neutron star matter will be rather small 
compared with the other isoscalar fields. 



In addition, the importance of off-shell corrections has been pointed out in detail in |33 
The scalar term has to change sign when going from on-shell u>k — m K to the off-shell point 
io k = 0. In the one-boson exchange picture this can be incorporated by adding the following 
terms 

C" K = D*KD^K ( f — + 1^1) + D *KTD»K f — . (35) 
\ m K m K J * m K 

The on-shell constraints eqs. (^) and (|33|) are still fulfilled if one replaces g„K by g„K — faK 

and gsx by g$K — fsx, respectively, for the values given in fl3"4|) . The off-shell behaviour of 

the scalar contribution is now proportional to 

I , faK u K \ g a K9uN 

oc 1 T — m K 

\ 9aKm K J mi 



and changes sign for uj = vn^j y/2, as it should be [S3], if one takes f a K — ^9aK- Hence, the 
decrease of the effective mass of the kaon in the nuclear medium is reduced by these new 
terms 

* 2 rn 2 K + g aK (rm K 
1 + faKcr/m K 

and makes kaon condensation even less favourable at higher densities (note that both the 
a-field and the coupling constants are negative here). 

We have found that the fit based on the adjustment to the KN-scattering lengths leads 
to an optical potential around U^ t = —(130 4- 150) MeV at normal nuclear density for the 
parameter sets used. This is between the two families of solutions found in . 
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The general expression for the effective energy of the charged kaon and the antikaon in 
neutron star matter is given by 

^k+,k- = m K J — ■ — jr ± {g uK Vo + 94>k4>o + 9 p kRo,o) (37) 

V m K + 2(g aK <7 + g a * K a*±g SK d) 

and shown in fig. 8 for the parameter sets TM1 and TM2 and model 2 including the 8- 
meson. The two cases (with and without the off-shell terms) are plotted. The energy of 
the kaon is first increasing due to the low density theorem. The energy of the antikaon 
is decreasing steadily. After hyperons are present the situation changes dramatically. The 
potential coming from the 0-field cancels the contribution coming from the w-meson. Hence, 
at a certain density the energies of the kaons and antikaons are equal to the effective mass of 
the kaon, i.e. the curves for kaons and antikaons are crossing at sufncent high densities and 
the energy of the kaon gets even lower than that of the antikaon! The difference between 
the two parameter sets considered is not as important as the inclusion of the off-shell terms. 
These terms effectively shifts the effective mass of the kaon to higher values. Therefore 
the in-medium energy of both kaons and anitkaons is shifted up to at most f 00 MeV. The 
minimum energy of the antikaon reads about 270 MeV for the on-shell and 350 MeV for 
the off-shell case. As the chemical potential never reaches values above 160 MeV here (or 
200 MeV without the 5-meson) antikaon condensation does never occur. We have checked 
the possibility of antikaon condensation for all parameter sets for model 1 and model 2 and 
found that at least 100 MeV are missing for the onset of condensation. The most uncertain 
term for the effective energy of the antikaon is the one coming from the <r*-exchange. If 
one arbitrarily doubles the coupling constant this would lead at maximum to an additional 
attraction of about 70 — 80 MeV and still no antikaon condensation is possible within our 
approach. 



V. CONCLUSION AND OUTLOOK 



We have studied the equation of state of neutron star matter within the relativistic mean 
field model including hyperons. We have used a representative set of parameters fitted to 
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the properties of finite nuclei and included the rather strong hyperon-hyperon interaction 
induced by additional (hidden) strange meson fields. The coupling constants of the hyper- 
ons have been fixed by SU(6)-symmetry relations and hypernuclear properties. Standard 
parametrizations with different forms of the scalar selfinteraction yield negative values for 
the nucleon effective mass solely due to the abundant presence of hyperons at high densities. 
The vector selfinteraction terms and the additional meson exchanges together eliminate this 
behaviour. According to our calculations, the baryon effective masses become small in neu- 
tron star matter at very high densities (p > 6po) and the abundances of the different species 
are determined by their spin-degeneracy factors (i.e. isospin-saturated matter). In all cases 
the hyperons get into the game at quite low densities, around 2 — 3p , and the composition is 
quite similar for all parameter sets studied up to 4p . The effects coming from the inclusion 
of a scalar isovector meson a (980) have been discussed. The equation of state does not 
change considerably but all the baryons get now different effective masses. The protons 
and electrons are getting a little bit suppressed at lower densities and the elctrochemical 
potential decreases while the S^'s occur at slightly lower densitites compared to the model 
without the 5-meson contribution. 

The possibility of antikaon condensation has been also studied in the framework of the 
relativistic mean-field model. Fixing the vector coupling constants by SU(3)-constraints, 
the in-medium energy of the antikaon was calculated for different choices of the optical 
potential at normal nuclear density. The existence of hyperons, especially the repulsive 
contribution from ^-exchange and the decrease of the electrochemical potential, makes a 
antikaon condensed phase less favourable even at very high densities. We get different 
answers for different parameter sets whether or not antikaons are present in the dense interior 
of neutron stars. Nevertheless, we see no onset of antikaon condensation below p < 6po for 
all cases studied. Note that the maximum density of neutron star will set another limit for 
the kaon condensed phase which has not been studied here. 

In a second approach we fix the coupling constants to the available KN-scattering data 
and take care of the off-shell behaviour of the kaon energy. Hence, this approach fulfills the 
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low-density theorem for kaons. The effective energy of the kaons and antikaons are dictated 
by the behaviour of the effective mass of the kaon as the contributions from the vector fields 
cancel each other. For very high densities (p > 7po) if is possible that the effective energy 
of the kaon is even lower than that of the antikaon as the contribution form the (hidden) 
strange meson field exceds the one coming from the cu-meson field. The effective energy of 
the antikaon turns upward again at this density. Hence, antikaon condensation can never 
occur in any of our parametrizations as the electrochemical potential is always lower than 
the minimum energy of the antikaon. 

The extrapolation to hyperon-rich matter is associated with large uncertainties due to 
the unknown hyperon-kaon interaction (e.g. A'H-interaction). The effective energy of the 
kaon is very much influenced by large cancellation effects coming from the vector fields. A 
slight increase of the coupling constant might therefore change the antikaon energy by 100 
MeV which are at least missing for the onset of kaon condensation. Hence, there might be 
still some niches for antikaon condensation. Nevertheless, the overall trend is quite clear: 
the presence of hyperons makes the onset of antikaon condensation quite unlikely to happen. 

The vector selfinteraction terms as well as the hyperon-hyperon interactions might in- 
fluence also the maximum mass, the rotational frequency and cooling properties of neutron 
stars. These questions will be studied in a forthcoming work. 
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TABLES 

TABLE I. The coupling constants of the parameter sets used. The vector coupling constants 
for the hyperons are be taken from SU(6)-relations. The nucleons do not couple to the a*- and the 
0-meson (g a *N = g<f>N = 0). The A does not couple to the isovector fields (g p \ = 0). The coupling 
constants for the E's are the same as for the A except for the isovector coupling constant which 
is 5pS = 2<? p 7v. The parameters for the scalar and vector selfinteraction terms are not given, they 
can be found in the corresponding references. 
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9aN 


10.0553 
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10.4262 
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11.4694 


9wN 


12.9086 


12.9450 


13.3415 


12.8861 
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14.6377 


9pN 
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9aA 


6.23 
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3.45 


3.59 


3.52 


3.43 


3.49 


3.94 


9ujE 


4.30 


4.31 


4.45 


4.30 


4.20 


4.88 


9p3 


4.85 


4.38 


4.56 


4.81 


4.63 


4.68 
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9(f>S 


-12.17 


-12.20 
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TABLE II. The coupling constants for the kaons. The coupling constants to the a- and <5-meson 
are fixed by the s-wave KN-scattering lenghts. The vector coupling constants are chosen from 
SU(3)-relations. The coupling constant to the cr*-meson is taken from /o-decay [31]. 



Set 


nl-z 


nl-sh 


pl-z 


pl-40 


tml 


tm2 


Ref. 


14] [ 


15 


1 1 § 1 1 



g uK 1.85 2.05 2.20 2.27 1.93 2.27 

g uK 3.02 3.02 3.02 3.02 3.02 3.02 

g pK 3.02 3.02 3.02 3.02 3.02 3.02 

g a * K 2.65 2.65 2.65 2.65 2.65 2.65 

g^K 4.27 4.27 4.27 4.27 4.27 4.27 

g SK 6.37 5.59 5.89 6.31 5.87 5.94 
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FIGURES 



FIG. 1. The composition of neutron star matter with hyperons in model 2 using param- 
eter set PL-Z. The jump in the curves is due to the negative effective mass of the nucleon. 

FIG. 2. The composition of neutron star matter with hyperons in model 2 using param- 
eter set TM1 with vector selfinteraction terms. The matter approaches isospin-saturation 
at high densities. 

FIG. 3. The effective masses of the baryons versus the density for model 2 using set 
TM1. All effective masses remain positive for the density range considered. 

FIG. 4. The field potentials and the electrochemical potential versus the density for 
model 2 using set TM1. Quite high values are reached indicating large cancelation effects 
at high densities. 

FIG. 5. The equation of state in model 2 for various parameter sets. The causal limit 
p = e is also shown. 

FIG. 6. The effective masses of the baryons versus the density for model 2 using set TM1 
with 5-mesons. All baryons get a different effective mass due to the additional meson. 

FIG. 7. The effective energy minus the chemical potential of the K~ over the density 
for model 2 using various parameter sets. The upper curves are calculated for an optical 
potential of U*^ = -120 MeV, the lower ones for [/* = -200 MeV. 

FIG. 8. The effective energy of the kaon and the antikaon for model 2 with 5-mesons 
using set TM1. The electrochemical potential is also plotted. Antikaon condensation does 
not occur over the whole density region considered. 
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Neutron star matter in the RMF model 
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